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Quantum gravitational effects in loop quantum cosmology lead to a resolution of the initial sin- 
gularity and have the potential to solve the horizon problem and generate a quasi scale-invariant 
spectrum of density fluctuations. We consider loop modifications to the behavior of the inverse scale 
factor below a critical scale in closed models and assume a purely thermal origin for the fluctuations. 
We show that the no-go results for scale invariance in classical thermal models can be evaded even 
if we just consider modifications to the background (zeroth order) gravitational dynamics. Since a 
complete and systematic treatment of the perturbed Einstein equations in loop cosmology is still 
lacking, we simply parameterize their expected modifications. These change quantitatively, but not 
qualitatively, our conclusions. We thus urge the community to more fully work out this complex 
aspect of loop cosmology, since the full picture would not only fix the free parameters of the theory, 
but also provide a model for a non-inflationary, thermal origin for the structures of the Universe. 



I. INTRODUCTION 

The possibility that primordial thermal fluctuations 
might seed the structure of our Universe is an intrigu- 
ing alternative to quantum fluctuations in a deSitter 
phase [l], H, H, 13] • Unfortunately a number of obstacles 
present themselves to such an enterprise. Firstly any 
thermal scenario should necessarily be based on a solu- 
tion of the horizon problem. This is so that the assump- 
tion of thermalization itself makes sense: modes must 
start off in causal contact to thermalize and then leave 
the horizon. (This is in fact true of any scenario where 
the fluctuations are "passive" @, @.) A number of solu- 
tions to the horizon problem have been proposed 0, ll> 
and in this paper we use in effect a combination of two. 

But more importantly thermal scenarios run against 
an apparent wall: the well know fact that thermal fluc- 
tuations have a white-noise spectrum, i.e. spectral index 
ns = 0, rather than the observed near-scale-invariance, 
ns ~ 1. Thus any scenario where the primordial fluc- 
tuations result from a "snap-shot" of a thermal bath at 
a fixed temperature is doomed. This discouraging re- 
sult, however, may be circumvented by noting that the 
white-noise nature of thermal fluctuations follows from 
the extensive nature of the energy. Reasonable and gen- 
eral as this assumption might be, it could be violated in 
the early Universe, during a phase ruled by new physics 
at the Planck or string energy scale. This has been sug- 
gested by at least two lines of research. In one a gas of 
strings at the Hagedorn phase is employed 0. Another 
invokes a holographic phase in loop quantum gravity [l(J . 

Yet another solution is to ensure that different modes 
leave the horizon and freeze-out at different tempera- 
tures. Then, the equal-temperature spectrum of ther- 
mal fluctuations is still white-noise, but the spectrum of 
frozen-in fluctuations imprinted outside the horizon isn't. 
The actual form of that spectrum depends on the bal- 
ance between the size of the mode leaving the horizon at 



a given time and the temperature (and thus the mode's 
amplitude) at the time the mode is picked out of the 
thermal bath, leaves the horizon and freezes-in. There is 
some controversy over whether this mechanism may lead 
to a scale-invariant spectrum and one of the purposes of 
this paper is to clarify the matter. 

In Section |TT] we provide a model calculation based on 
a minimally modified thermal scenario, in which thermal 
matter is allowed to have a different equation of state 
with w = p/p < —1/3, but where nothing else is changed. 
We show that unless new physics comes into play, mod- 
ifying the Einstein equations, the thermodynamical rela- 
tions, or some other standard assumption, in all such 
scenarios the spectral index is ns — 4. This is true re- 
gardless of w, the only free parameter of the model. This 
section is supplemented by Appendix [X] which provides 
all relevant definitions of measures of structure. 

Of course it is natural that new physics does come 
into play in the early Universe, and the rest of this paper 
is focused on the potential of loop quantum cosmology 
(LQC) to reverse this negative result (for an up to date 
introductory review see and for early developments 
in the field see [13]). Modifications to the the Einstein 
equations in LQC originate from two sources: the field 
strength (curvature) of the Ashtekar connection which is 
expressed in terms of holonomies and the inverse pow- 
ers of volume in the constraint. These are quantified 
by two parameters in the theory. Firstly, the j param- 
eter, which appears due to tracing over holonomies of 
the SU(2) connection in both gravitational and matter 
parts of the constraint. This parameter also determines 
the scale at which modifications to inverse volume be- 
come significant. Secondly, the I parameter, which arises 
due to the inverse volume term in the matter part and 
quantifies the functional form of this modification. 

So far a complete and consistent quantization of LQC 
with the knowledge of physical Hilbert space has only 
been performed for j — 1/2. It shows a physical resolu- 
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tion of the singularit y a nd leads to the correct classical 
limit [IE EE EE EE 111- Nevertheless, these investiga- 
tions have provided valuable insights on existence of ef- 
fective Hamiltonian which is an excellent approximation 
to the underlying quantum dynamics and which can be 
generalized to higher j in particular for the regime where 
the Hubble rate is small compared to Planck scale. In 
this regime, modifications to the Einstein equations are 
primarily due to corrections to the inverse volume opera- 
tor. Though these modifications unfortunately cannot be 
tied to a curvature/ energy density scale in non-compact 
flat models [IE E3 > they are related to the scale of in- 
trinsic curvature in closed models and as we will show 
also lead to a bounce of the closed universe at a critical 
scale factor [H[. Our analysis will be carried out in this 
framework. 

In Sec IIIII we present this effective Hamiltonian and 
the derivation of the background dynamical equations. 
To keep the paper suitable for a phenomenology oriented 
readership, all the necessary details of loop quantization 
are relegated to Appendix [Bl where we give particular 
emphasis to the origin of the j and I parameters with 
resulting correction terms in the new improved quanti- 
zation of LQC 0, EE El, El and the role of differ- 
ent modifications to the Einstein's equations. Equipped 
with the background dynamics in LQC, in Section HVl we 
then investigate the spectrum of thermal fluctuations. 
We demonstrate that the modified dynamics in LQC 
evades the no-go result of Sec. HU The constraint on 
the near scale invariance of fluctuations translates into 
the requirement on the effective equation of state in loop 
cosmology, w « —2/3. This can be obtained by a suitable 
choice of the I parameter. 

Since work on inclusion of inhomogenities in loop cos- 
mology is at a very early stage of development [HI (and 
various technical aspects, in particular those relevant for 
the regime of interest, are yet to be understood) we adopt 
a phcnomcnological approach and parameterize the ex- 
pected modifications to the perturbed Einstein equa- 
tions. We show that expected corrections in fact make 
the above result stronger. 

More work needs to be carried out by the community 
until our calculations may be converted into specific con- 
straints. But, as we summarize in SectionfVl in this paper 
we are able to provide a list of what exactly needs to be 
worked out so that n$ ~ 1 is converted into a constraint 
upon the free parameters of LQC. 



II. THE BASE CALCULATION AND THE 
NO-GO RESULT IN CLASSICAL PHYSICS 

Let us consider a thermal scenario in which the only 
effect of new physics is to change the equation of state of 
thermal matter. This certainly happens in theories with 
deformed dispersion relations and also in LQC [20| . 

We then assume that statistical physics, the gravitational 
equations and a few basic thermodynamical relations are 



not modified. This won't necessarily happen in loop cos- 
mology, but the way in which the modifications arise is 
not yet fully understood. It is therefore interesting to 
provide a base calculation, assuming no chang es, as a 
blueprint for further work. 

The calculation follows three steps. 



A. The fixed temperature power spectrum 

The first step is to compute the fixed temperature 
power spectrum. This turns out to depend only on the 
specific heat at constant volume (a result that has been 
known since the XIX century, and is now textbook ma- 
terial [21]). The spectrum is generally white noise, a fact 
that can be directly traced to the extensive nature of the 
energy, i.e. to the fact that the energy inside a given 
region is proportional to its volume. 

The derivation is very general. Consider the partition 
function 



Z = J2 ( 
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(1) 



where (3 = T 1 . The total (matter) energy U inside a 
volume V is given by: 



U=(E) 
and its variance by 
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where cy is the specific heat at constant volume. If 
the energy is extensive then U = pV, with energy den- 
sity p = p(T), that is U is proportional to the volume. 
The spectral index and amplitude may now be found by 
means of the following tools (for the closed model we 
will be interested in the regime where intrinsic curvature 
terms can be neglected): 



• The Poisson equation 

/c 2 $ = AnGa 2 p8 = AnGa 2 5p, 



(4) 



relating the gravitational potential and the den- 
sity fluctuations. Outside the horizon there may 
be gauge issues, but this relation certainly holds 
for subhorizon modes. 

• The proportionality between the variance a 2 (R) in 
a quantity g defined in position space and smeared 
on a scale R, and the "dimensionless power spec- 
trum" : 



a: 



:{R) = (Sg 2 



V g (k R = a/R) 



(5) 



(see Appendix [Al for definitions; note that here k 
is comoving, but R is a proper size). The spectral 
index ns is defined from V$ — A 2 fc" s_1 . Formula 
([5]) has not been questioned for ns < 1; but see [22l . 
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Then using Q we have 



Combining this result with ([4]) and 
thus conclude 



- R ecv ■ (6) 
(for g = Sp) we 
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F 



Cy 



= ^ T2 M«=f ■ (7) 



Using the extensive nature of the energy we have cy = 
p'{T)R 3 , so finally 



-T 2 p' . 



(8) 



The fluctuations are therefore white noise (n<j = 0), and 
have an amplitude that only depends on the Stephan- 
Boltzmann law, relating energy density and temperature. 
We shall assume a Stephan-Boltzmann law of the type 
p oc P, where £ is a parameter subject to thcrmody- 
namical constraints to be discussed later. 



B. Frozen-in power spectrum 

We assume that comoving scales k start thermalized 
and inside the horizon, and then leave the horizon, with 
(first) crossing defined by k = aH. This requires either 
accelerated expansion Q , a loitering phase [|[ , a decreas- 
ing speed of light 0], a bouncing scenario, or a combi- 
nation thereof. We use the first mechanism, so that the 
equation of state satisfies w < —1/3. 

As the <!> modes leave the horizon their amplitude gets 
fixed at whatever thermal amplitude they have at cross- 
ing, that is: 



-T 2 p' 



k=a.H 



(9) 



Since different modes freeze at different temperatures the 
spectrum left outside the horizon won't be white noise. 

Using the Friedman equation H 2 oc p we can rewrite 
© as 



T 2 p' 



(10) 



k=aH 



where k = aH specifies a relation between a given co- 
moving k leaving the horizon at a given time, and the 
temperature, thereby allowing the inversion of the right 
hand side as a function of k. Eqn. (fTO)) implies: 



= 1 + - 

dlnT 2 



(11) 



The relation between k — aH and the temperature, how- 
ever, depends on both the equation of state p = wp 
and p oc Using the Friedmann equation we have 

k = aH oc a^fp, and since p oc l/a 3 ^ 1+w \ we may derive 



a oc Tv+») . 



(12) 



Therefore: 

dink -C C _ C(l + 3w) 

dlnT ~ 3(1 + w) + 2 ~~ 6(1 + w) 

We can now compute the spectral index as 



(13) 



dlnT® dlnV^dlnT 2 + (l + w , ^ 
ns-1 = , = ^ , = 3— — - ; - (14) 



din A: 



d In T d In k 



( l + 3w 



but note that the condition w < —1/3 (or that k = aH 
increases in time) is necessary for this formula to make 
sense. 

Two promising regions of parameter space stand out. 
Firstly £ = — 2, that is p oc 1/T 2 ; this may lead to scale- 
invariance because the amplitude of the frozen-in thermal 
fluctuations does not depend on the temperature in this 
case (c.f. Ean fTO)) . Secondly w = — 1; one can see that 
this could lead to scale-invariance because p does not 
change (it behaves like a cosmological constant), and so 
neither does the temperature or amplitude of the fluctu- 
ations as they leave the horizon. 

However further conditions apply. Regarding the first 
case we have to check that w < — 1/3 is possible, so that 
modes do leave the horizon. With respect to the latter, 
we should additionally have £ ^ (or £ ^ oo), so that 
there are fluctuations at all (and they are not infinite). 
Unfortunately closer inspection shows that these condi- 
tions cannot be met. 



C. Thermodynamical constraints 

It's been noted [H, [2f| [2(| that the equation of state 
p = wp and the Stephan-Boltzmann law p = p(T) are 
linked by a thermodynamical relation. The argument 
assumes that energy and entropy are extensive. Consider 
the first law of thermodynamics: 



dU 



-PdV + TdS 



(15) 



If the energy U and entropy S are extensive, then 
U(XV, XS) = XU(V, S). Taking a derivative with respect 
to A at A = 1, and using (fl5| ) we arrive at the Euler rela- 
tion 



U = -PV + TS 



(16) 



so that defining p = U/V and entropy density s = S/V 
we have 



P + p 



(17) 



We can now prove that s = dP/dT in a variety of ways, 
e.g. introducing the free energy F = U — TS = F(V, T), 
so that dF = -PdV - SdT. This leads to the integra- 
bility condition: 



dS_ 

dV 



dP 
df 



(18) 
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Thus the expression 



dP 
dT 



T 



If w is a constant we obtain that p oc T< with 



(19) 



(20) 



The trouble is that this relation implies that £ = for 
w = — 1: "deformed" radiation may behave like a cos- 
mological constant, but then the specific heat vanishes 
and there are no thermal fluctuations at all. This is an 
interesting result but kills the second candidate for scale- 
invariance proposed above. 

The first candidate is killed by noting that ( — —2 im- 
plies w = —1/3, that is a Milne Universe. This is merely 
a borderline case for solving to the horizon problem: the 
comoving horizon does not increase but neither does it 
increase. 

One might expect that models near these two can 
bypass these problems and display if not strict scale- 
invariance, at least approximate scale-invariance. How- 
ever this is not the case. Inserting (|20]l into (fT4|) a simple 
algebraic calculation shows that w (or £) cancel out and 
we are left with 



n s 



(21) 



for all model parameters. 

Therefore one needs further new physics to bypass this 
negative result. Presumably the double branched disper- 
sion relations considered in 0] are behind the fact that 
w = — 1, C = 1 is possible (in contradiction with (|20| . 
We now examine the way modifications to the dynamics 
in LQC reverse these results. 



III. EFFECTIVE DYNAMICS IN LOOP 
QUANTUM COSMOLOGY 

The phase space in LQC consists of the geometrical 
variables - the connection c and the triad p - and the 
matter variables, which for a scalar field will be <fr and its 
canonical momenta p$. The triad is related to the scale 
factor as p = a 2 = V 2 ! 3 . On the classical solutions of 
GR for closed model, c is related to the time derivative 
of scale factor as c = ja + 1 where 7 w 0.2375 is the 
Barbero-Immirzi parameter. The connection and triad 
are canonically conjugate satisfying 



{c,p} = 



8nGj 



(22) 



For the closed model, in the regime where the Hubble 
rate is small compared to Planck scale we can write an 
effective Hamiltonian which encodes the modifications to 
the inverse scale factor below a critical scale a* (param- 
eterized by j (|B2ip ) in terms of functions S (Ea. (|B22p ) 
and Di (Eq.jB23|): 



H c ff — — - 



8irGr 



■Sa ((c- 1) : 



'7 



(23) 



where H m is the matter Hamiltonian obtained after in- 
verse volume modifications using (|B20|) . For a massive 
scalar field it is 



2 

H m = \dM + a 3 V{<j>) . 



(24) 



Dynamics can now be obtained by the use of Hamil- 
ton's equations. In order to obtain the modified Fried- 
man equation we first evaluate 



P = 



{p, H cS } = - 
2Sa , 

and then using it in H c g ps 0: 



SttG H„ 
8ttG 



"~3 IhT 



(25) 



9 

a 
a 2 



(26) 



where p. Ha = H m /a 3 denotes the modified energy den- 
sity [23] (we follow the conventions of Ref. [13]). It can 
be easily seen that a bounce occurs for a < a* when 
(87rG/3)p sc = S/a 2 which is possible due to form of S 
and Di in this regime. Also, for a » a*, 5 and D\ ap- 
proach unity yielding us the classical Friedman dynamics. 

We can obtain the modified Raychaudhuri equation 
using Hamilton's equation for c: 



{c, H eS } 



8ttG7 dH cS 



3 dp 

and the expression for the modified pressure 
dH M 



p 

dV 

These equations lead to 



2 _y 2 dH M 

3 dp 



(27) 



(28) 



4ttG7 aP sc (29) 



~ ((c-l) 2 +7 2 ) - + - 
27 \ a a 



which combined with the time derivative of Eq. (|25p result 
in the modified Raychaudhuri equation: 



a 4ttG „, „„ . IS fa S 

- = —S( Psc + 3P sc ) + -- 

a 6 I a \o a 



(30) 



It is then straightforward to verify, using the Friedman 
and the Raychaudhuri equations, that p satisfies the con- 
servation law: 



Psc + 3-(Psc + Psc) = 



(31) 
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Defining the modified equation of state as 

W sc = Psc/Psc (32) 

we therefore have 



Ps 



POs 



a 3(l+w. 



(33) 



The modified Klein-Gordon equation can be derived 
by using 



= {<t>,H m } = 



dH m 

dcj> 



and 



, H m } — — 



dp4, 



(34) 



(35) 



Taking the time derivative of (p we are led to 



d / 1 din A 
a V 3 d In a 



dV 

D ' in ■ ( 36 ) 



For a <C a*, Di <C 1 and the dynamics of the massive 
scalar field just behaves as of the massless scalar field. It 
is straightforward to verify that the energy density 



^ 4-V 



and pressure 



SC "2D, I 3 dlna ' 



(37) 



(38) 



lead to above Klein-Gordon equation through the con- 
servation equation (j3Tj) . 

The modified equation of state thus becomes 



P,, 



(39) 



Psc (f> 2 + 2Di V 

Using Eq. (|B23|) we obtain w sc for a <C a*: 

w sc ^l-a (40) 

where 

3-Z 



1 - Z 



(41) 



The modified equation of state for arbitrary matter can 
be similarly found by following the procedure in Ref. [2(| • 
Here one views p sc as being obtained from substituting 
inverse powers of scale factor as appropriate powers of 
Di . If the classical energy density is given by 



Po 



then it is easy to see that 



a 3(l+w c ) 



D" 



(42) 



(43) 



Since the latter satisfies a conservation equation, it 
evolves according to an expression like (|4"2"|) but with 
modified equation of state 



w sc = w c 1 



1 din A 
3 din a 

Thus for a <C a* we obtain 

w sc « W c (1 - a) . 



(44) 



(45) 



Since < I < 1, w sc can be easily less than —1/3 for 
arbitrary matter when o < a*. 



IV. ORIGIN OF THERMAL FLUCTUATIONS 
IN LOOP COSMOLOGY 

In the preceding section we saw that even at the zeroth 
order there are modifications to the Friedman dynamics 
in loop cosmology. These are sufficient to possibly over- 
come the no-go result obtained in classical physics. We 
start with the simplest possibility where the only change 
from the classical physics appears via (f2"fj|) and (|3H)) . As 
in Sec. II, we will be interested in the scales where the 
modifications coming from the intrinsic curvature can be 
ignored in the Friedman equation. We note that immedi- 
ately after the bounce, p sc becomes dominant over S/a 2 
term and for a proper choice of initial conditions for mat- 
ter such a regime can coexist with a <C a*. 

Using ©, Ea.([T0]) modifies to 



2 „/ 1 



T 2 p 



VSpl 



Since p sc oc a 3 ' 1+B ' C ' and p cx we obtain 



a oc T3(i+«.e) 



(46) 



(47) 



Using Eq. (|B22[ ) we find that for a -C a*, S cx a 3 and 
hence 



dlnP* (l + w sc )(2 + C) + C 



d\nT 



2(1 + w sc ) 



(48) 



Also fc = aH oc a^/Sp sc leads to 

dink C(l + 3w sc )-3C 



dlnT 6(1 + w sc ) 

The modifications to the spectral index thus become 
dlnV^ „ (l + w sc )(2 + C) + C 



(49) 



n s - 1 



dlnfc C(l + 3w sc )-3C 



(50) 



In addition if we use a semi-classical density for the en- 
tropy, 



S_ _ £[_ 
V ~ ^ 



(51) 
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the thermodynamical argument presented in Section [TH 
relating £ and w, is also valid, for the semi-classical values 
of these parameters. Thus, 



Wsc 



(52) 



On substituting this in Eq. (|50|) . the condition for near 
scale invariance ns ~ 1 translates to the requirement 
that w sc ~ —2/3 which can be obtained by an appro- 
priate choice of I parameter and w c 28]. Therefore we 
find that at the zeroth order the no-go conditions for 
scale invariance of thermal fluctuations can be overcome 
by modifications to the gravitational dynamics in loop 
cosmology. 

However, the calculation presented above is incom- 
plete, because very little is known about the perturbed 
Einstein's equations in loop cosmology, in particular for 
a « a t . Preliminary work on flat models and the 
regime a 3> a* suggests that modifications to the gravi- 
tational dynamics influence the growth of fluctuations in 
a very non-trivial way. Based on these calculations we 
classify below some possible modifications: 

• In a simplified setting the Poisson equation inside 
the horizon could become 



fc 2 $ = 4nGa 2 D I (a)5p, 
where I is an unknown exponent. 



(53) 



• It could be that the scale where the fluctuations 
become dominated by gravity (and not pressure) is 
not simply proportional to the horizon scale k ~ 
aH . For simplicity we shall ignore this possibility: 
it relates to varying speed of sound scenarios to be 
explored elsewhere. 

• It could be that beyond the gravity-driven "freeze- 
out" scale the potential $ continues to evolve and 
does not freeze-out as usual. This was proved ex- 
plicitly in [L9J for a > a*. Here we shall model the 
evolution of the potential outside the horizon for 
a < a+ as 



N 



$ oc a 



where N is an exponent to be computed. 



(54) 



We stress that these modifications parameterize our igno- 
rance of the theory but they should be derivable in terms 
of j and I alone, from first principles. 

Given these novelties we find that formula |(9]) gets 
modified to 



tD 2I t 2 p ' s 

k 



k=aH 



and since H cx \JSp sc we have at horizon crossing 



T 2 D 2I p' sc 



(55) 



(56) 



But because the potential continues to evolve outside the 
horizon this is not enough to read off a condition for scale 
invariance. Indeed the spectrum left after a = a+ will be 
processed into 



T 2 D 2I p' sc 



VSps 



,27V 



(57) 



H 



where all the quantities on the right-hand side are to be 
computed when the mode left the horizon, for k = aH. 
The relation (j 1 2|) between temperature and a is unmodi- 
fied (apart from replacing £ and w by their semi-classical 
values). Therefore we can deduce the counterpart of (fTT|) 
as 



fgg» =1+ Cf 1+ 1 

d\nT H 2\ 3{l + w sc ) 



(3 - 12/a + 4iV) 



where Th is the temperature when the mode left the 
horizon. 

Using Eq. (|4"9")) we then obtain 

3(l + Wse )(C + 2) + C(3-12/ Q + 4iV) 
C(l + 3w sc ) - 3C 

The condition for near scale invariance (ns ~ 1) then 
implies 3(1 + w sc )(( + 2) w -((3 - 127q + 47V) leading 
to 



2 4 

3 ~ 9 



(N - 3/a) 



(60) 



iH 



This equation can be viewed as a first order improve- 
ment over our zeroth order calculation which led to 
w sc w -2/3. 



V. CONCLUSIONS 

Loop quantum cosmology has the potential to relate 
observational physics and quantum gravity, allowing con- 
crete calculations to be made in the quantum gravity 
regime as long as a minisuperspace approximation is as- 
sumed to be valid. The approach is known to modify the 
equation of state of ordinary matter, thereby permitting 
a solution of the horizon problem without resorting to 
scalar fields. It is then natural to ask whether in such sce- 
narios thermal fluctuations could be behind the observed 
structure of the Universe. In order to analyze this issue 
we have assumed in this paper that physics learned in the 
mini-superspace approximation (in the sense of modifica- 
tions to inverse volume terms) will not change qualita- 
tively. There are positive indications for this hope from 
ongoing work (29| , but we stress this important caveat in 
our analysis. 

We showed that prima facie we are confronted by a 
no-go result in classical physics, pointing to ns = 4 in 
all such scenarios. This can be derived assuming only 
the Einstein equations and a basic thermodynamics re- 
lation. The fact that the zeroth order, background, Ein- 
stein equations are also modified in loop cosmology allows 
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us to bypass this negative result, pointing to the region 
of parameter space where ng w 1 is realized. This oc- 
curs for the semi-classical equation of state w sc « —2/3. 
However, before this requirement can be converted into 
a constraint upon the free parameters of the theory (J 
and I), a number of important details have to be worked 
out. We close with an executive summary of what is still 
missing in the theory: 

• A solid quantization in the regime of large j, neces- 
sary for a full understanding of an extended period 
with a < a*. This has to be accomplished for closed 
models as for non-compact flat models the physical 
meaning of a* makes little sense. 

• A study of the perturbed Einstein equations along 
the lines of that carried out in [19(, but valid for 
a < a*. 

• A concrete prediction for the spectrum of gravita- 
tional waves (tensor modes) completely ignored in 
this paper. 

We believe that once this task list is completed we shall 
be able to place solid observational constraints upon loop 
quantum cosmology. 

We conclude with a final remark on the role of higher 
j terms in both the gravitational and the matter parts 
of the Hamiltonian. In various LQC phenomenology pa- 
pers one has often ignored the modification to the gravi- 
tational part (constituted by S) . Such an ad-hoc analysis 
is similar to taking different metrics in the gravity and 
matter parts of the Einstein equations in GR. As a purely 
academic exercise we can perform such an analysis and 
it turns out that the zeroth order calculation in Section 
[TV]does not go through; instead one reproduces the no-go 
result for scale invariance of classical physics [3lJ . But by 
consistently incorporating modifications arising for high 
values of j in both the gravitational and matter parts of 
the constraint, the no-go obstacle is removed even at the 
zeroth level of calculation. This is an important lesson for 
loop cosmology phenomenology, showing the non-trivial 
features of high j. We believe this opens an interesting 
avenue for re-examining various interesting ideas (e.g., 

Refs. PUIS). 



APPENDIX A: MEASURES OF STRUCTURE 

For discrete Fourier modes we define, for any quantity 
g, the dimensionless (or curly) power spectrum V g (k) as: 



Vg{k) 



V 



;47rfc 3 (| 5k | 



BK "' (2tt)3 
or occasionally the non-curly one as: 

P 9 (k) = V{\g k \ 2 ) = ^-V g (k) 



(Al) 



(A2) 



The latter is often only used for 5 and using the Poisson 
equation we have the following alternative definition of 
the the spectral index 



Ps(k) — A 2 k r ' 



(A3) 



The Fourier transform can be introduced noting that 
Afc = 2tt/L so that: 



v 4^ 



Then with: 



S(x) 



(2tt)3/2 



dk g(k)e 



zk-x 



6{k) = (2^/ dX ' 



S(k) = 



(2tt) 3 / 2 



we have 



5k ~ s(k) 



dx g(x)e~ 



(2tt) 3 / 2 



V 



<W ~ y 0{k - k ) 



(A4) 



(A5) 
(A6) 
(A7) 

(A8) 
(A9) 



and so we find the alternative and equivalent definition 
for the power spectrum: 



2tt 2 

(g(k)g*(k')) = —V g (k)8(k-k') 



(A10) 



The position space variance, with either definition, can 
be written: 



Acknowledgments 

We thank Stephon Alexander for useful discussions in 
early period of this work and Lee Smolin for helpful com- 
ments. Research at PI is supported in part by the Gov- 
ernment of Canada through NSERC and by the Province 
of Ontario through MEDT. PS is supported by the NSF 
grant PHY-0456913 and the Eberly research funds of 
Penn State. 



^ = <^(x))= / ^r g (k) 



(All) 



The filtered position-space variance is also used. It's 
based on the smoothed field 



g(R,x) = — j g{x')W{\x-x'\/R)d^' (A12) 



V R = / d A xW{x/R) = AnR A / y z W{y)dy{M 3) 
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where R is the smoothing scale and W can be, say, a 
Gaussian or a top hat. (Do not confuse Vr with the 
large V used in the discrete Fourier series.) Then the 
"sigma-squared" on scale R is 



a 2 JR) = (g 2 (R,x)) = 



dk 



W 2 (kR)V g {k) (A14) 



where we used the convolution theorem and W(kR) is 
the Fourier transform of (2tt) 3 / 2 W(x/R)/Vr, that is, it's 
normalized so that W(kR) = 1 at k = 0, then falling off 
at fcfl ~ 1 / R. We can then write approximately 

a*(R)ttV g (k R ) (A15) 

for kit = 1/L, since that's where the integrand peaks. 



APPENDIX B: SOME BASICS OF LOOP 
QUANTUM COSMOLOGY 

This appendix aims to summarize some key features 
of loop quantization of cosmological models (for de- 
tails see for example Refs.flJ, [la]). We first demon- 
strate the way the classical constraint is cast in terms 
of elementary loop variables - holonomies of connection 
h e (A) = Vexp(J e A) and triads Ef. We also show the 
way j and I parameters and the resulting modifications 
in the improved quantization of LQC [HI to classical GR 
originate. 

For simplicity we start with quantization of flat 
isotropic and homogeneous FRW spacetime. In this set- 
ting the underlying symmetries lead to simplified connec- 
tion c and triad p : 



Al 



cVo 



and Ef = p V 



fa~ °p a 

"do L 4 



(Bl) 



where (°L0 l al °e1) are a set of orthonormal co-triads and 
triads compatible with the flat fiducial metric °q ab . V a 
is the volume of the cell (V) used to define a symplectic 
structure with respect to °q ab - The variables c and p are 
canonical conjugate satisfying Eq. ([22|) . 

The gravitational and matter parts of the constraint 
are given by 



Cgrav T 



c^xe tjk e- 1 E ai E^F! e 



ah 



and 



8ttG 



\p\3/2 



(B2) 



(B3) 



where for simplicity we consider a massless scalar field. 
The modulus sign arises because of two possible orien- 
tations of the triad, the choice of which has no physical 
consequences unless we choose spinor fields. F* b denotes 
the curvature of connection and 



e = y/\detE\ = 



l ^ Z? a T?b rpc 

e ijk E t E m E n 



1/2 



(B4) 



To write the constraint in terms of holonomies and triads, 
the following identities of classical phase space are very 
useful [H: 



1 



87rG 7 



{Af, e^ k e abc EfE^El} = 3 e^ l e abd E?E$ , (B5) 



{AlV} {Ai,V^} 



yn 



and 



1 -n 



{AlV} 



(B6) 



(B7) 



Eq.(|B5|) leads to 

e ijk e- l E ai E b iF* b 



8ttG7 



e- 1 e abc {Al,V<} F abi 
e abc {A^V} F abl . (B8) 



Here V 



\P\ 



3/2 



a 3 denotes volume of the cell with 



respect to the physical metric V = V ^ v /\detE\ (for sim- 
plicity we put V a = 1 from now on). 

We then express the connection in terms of the holon- 
omy by tracing over the holonomies in a j representation. 
For j = 1/2, using 

Tr( n r j ) = -^j(j + l)(2j+l)5i , (B9) 



it is straightforward to obtain 
e ijk e- x E ai E b i = ]T 



,abc o. ,k 



2ttGj\ 



Tr(hi X) {hi X) -\V}n) 



(BIO) 

Here h k is the holonomy of the connection c along the 
edge A°e£ 



hi ' = cos - 1 + 2 sin — r k 
K 2 2 



(Bll) 



where I is the identity matrix, r k = —ia k /2 and a k are 
the Pauli spin matrices. 

The curvature components can be obtained by consid- 
ering holonomies around a closed square loop □jj-: 



F* b = -2 lim Tr 

Ar n ^0 



1 l Zl 



(B12) 

Due to the inherent quantum nature of geometry in loop 
quantization the area A 2 |p| is shrunk to the minimum 
eigenvalue of the area operator A = (2\Z3nj)£p in LQG 
[Tj]. This leads to a constraint A 2 |p| = A. C grav can 
then be obtained by combining (|B10p and (|B12[) . 
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The matter part contains inverse powers of det E. To 
write them in terms of holonomies we use the identity 
|dete l Q | = | det E\ x l 2 and Eq.jB7| 



det el 



6 det£ 



= e abc e ijk {Al V^}{Al V^ 3 }{A k cl V 1 / 3 } . 

(B13) 

Expressing connection in terms of holonomies, in j = 1/2 
representation inverse scale factor becomes 



sgn(p) 

|p|l/2 



8ttG 7 A 



Tr £ r k h 

k 



(A)-^ y l/3| 



(B14) 

Note that in obtaining (|B13|) we could have multiplied by 
(det e^/Vdet E) m where m has a value such that we ob- 
tain po sitive po wers of the volume in the final expression 
of l/yj det E\. This happens when m > 1/2. We are 
thus led to the / ambiguity in the expressions for inverse 
triads in the matter part of the constraint and Eq. (|B 14|) 
generalizes to 



sgn(p) 
|p|i/a 



l/2(l-i) 



(B15) 

where Z = 1 - l/(2m) and < / < 1. 

With C grav and C m written in form of holonomies and 
triads, we quantize the theory and are led to a non- 
singular difference equation with uniform discretization 
in eigenvalues of volume operator [lij : 



V\v) = [3 3 \v\\v) 



with 



V 6 



A- 1 / 3 , K 



2y/2 



3V373 



(B16) 



(B17) 



Unlike the old quantization in LQC (where the differ- 
ence equation was of uniform discretization in eigen- 
values of triad [H, [34|), the evolution has the correct 
classical limit for arbitrary matter content and quan- 
tum gravitational effects set in when curvature becomes 
of the order Planck. Study of backward evolution of 
semi-classical states peaked at late times on trajecto- 
ries of a large classical universe shows a generic bounce 
when p = p cr it = 0.82ppi [l4j]. In this quantization (for 
j = 1/2) modifications originating from F* b terms dom- 
inate over those containing 1/vdet E in both gravita- 
tional and matter parts of the constraint. Investigations 
of closed models yield similar results [l(| • 

We now provide the expressions of the eigenvalues of 
operators corresponding to e~ 1 E at E b ' 1 and 1/ydet E for 
higher j , which can be derived in analogy with (|B10[) and 



(|BT4]) . using Eq.(|B9"l). We have 
9 /3 3 K 2 / 3 



87r4 7 j(i + l)(2j + l) 



H i/3 r\v-2r 



r=-j 



(B18) 



and 



27 P 21 K 2 / 3 



I6n£ 2 7 lj(j + l)(2j + l) 



\v\ 1/3 V r\v- 

1 2 j + l ) m 



. 2r |2//3 



(B19) 

For higher j, Eqs. (|Bl"8f and ([519)) can be approximated 

by 



Sj=S(q)a, and dji(q) = Di(q)a 3 (B20) 



with 



q:= (a/a,) 3 , a^ = (2jf 5 0, 



(B21) 



2{(q+l) 3 -\q-l\ 3 ) 



- 3g((g + l) a -sgn(?-l)|g-l| 



(B22) 



and 
DM) 



27 Iqf-f r 1 



8/ 



r / 2(i- 



+ 3) 2(1+3) 

- g-1 3 



2q 



21 + 3 



/ 2(1 + 3) 2(1 + 3) \ "1 

({q + 1)-^ - sgnfa -l)\q- 1|^) } 



2(1-0 

(B23) 



following the analysis for old quantization [35|, [3g|. In 
the regime when a <C a* 



5(9) 



3 ( a 



2 W 



and 

A(<?) 



g ^3/(2(1-0) f a \ 3(3-0/(1-0 



2Z + 3 



(B24) 



(B25) 



For a ^> a*, 5(g) w 1 and ~ 1- 

Extensive numerical simulations of backward evolution 
of semi-classical states at late times have shown that in 
LQC an effective Hamiltonian which provides an excel- 
lent approximation to the underlying quantum dynamics 
can be written. For j = 1/2, the effective Hamiltonian 
for flat [3| and closed model [l|| is 



H, 



CW 3s j sin 2 (Ac) 

77^77 8^5^ + ^'Y (B26) 



and 



3a 3 - (sin(Ac)sin(A(c-l)) + (l+ 7 2 )) , , p£ 

Hcff = 8^ + j ' l ~2 

(B27) 
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respectively. Here A = X(p) — (A/lpl) 1 / 2 . The sin(Ac) 
terms arise from field strength part and are responsible 
for p 2 modifications of the Friedman equation [3, EB] • 

For higher j a quantization procedure has been pro- 
posed which indicates resolution of singularity [36j ; how- 
ever a complete quantization is still lacking. Considering 
higher values of j in non-compact flat models has the 
problem of relating the scale at which modifications to 
l/Vdet E terms become important to any physical scale 
(The scale at which these modifications become impor- 
tant depends on the choice of fiducial cell, for details see 
Refs.fT4.fT7j'). In the closed models this scale is provided 



by the intrinsic curvature. If we consider higher j in 
closed models then, as for j = 1/2, modifications arise 
both from F* b and l/V dot E terms. However, in this 
case we can have a regime in which modifications com- 
ing from the latter dominate the former and still lead to 
a non-singular bounce. When sin(Ac) — > A 1 / 2 c/a and 
sin(A(c — 1)) — > A 1 / 2 (c — l)/a, the modifications com- 
ing from F* b can be considered small. In this regime 
H = a/a < A" 1 / 2 ~ M Pl and energy density is small. 
An effective Hamiltonian (|2"5|) can then be obtained by 
following the procedure outlined for higher j {36} . 
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